Quarkonium dissociation in a PNJL quark plasma* 
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We investigate the Mott effect for heavy quarkonia due to Debye screen- 
ing of the heavy quark potential in a plasma of massless quarks and anti- 
quarks. The influence of residual color correlation is investigated by cou- 
pling the light quark sector to a temporal gauge field driven by the Polyakov 
loop potential. This leads to an increase of the Mott dissociation tempera- 
tures for quarkonia states which stabilizes in particular the excited states, 
but has marginal effect on the ground states. The temperature dependence 
of binding energies suggests that the dissciation of the charmonium (bot- 
tomonium) ground state by thermal activation sets in at temperatures of 
200 MeV (250 MeV). 
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1. Introduction 

Since the suggestion of J/ ip suppression as a signal of quark-gluon plasma 
(QGP) formation by Matsui and Satz [I] in 1986 the problem of quarkonium 
dissociation in hot and dense strongly interacting matter has played a key 
role for QGP diagnostics in relativistic heavy-ion collision experiments. The 
original idea was that in a QGP the string tension of the confining potential 
vanishes and the residual one-gluon exchange interaction undergoes a Debye 
screening by the color charges of the plasma. When the temperature depen- 
dent Debye radius r£>{T) (the inverse of the Debye mass m£>{T)) becomes 
shorter than the Bohr radius of the charmonium ground state (J/ip) then 
the Mott effect [2] (bound state dissociation) occurs and the corresponding 

temperature is T^ tt . This simple idea grew up to a multifacetted research 
direction when not only in the first light ion - nucleus collisions at the CERN 
NA38 experiment, but also in proton - nucleus collisions at Fermilab J/ip 
suppression has been found so that there is not only a QGP but also a cold 
nuclear matter effect on charmonium production, see [3] for a recent review. 

If one wants to explore the question of screening in a plasma more in 
detail then a variety of approaches is available in the literature, from the 
original Debye-Hiickel approach [5J applicable to any vacuum potential (for 
example the Cornell potential), over the thermodynamic Green functions 
approach to the ab-initio studies of heavy-quark potentials in lattice QCD. 
With the obtained medium-dependent potentials one can then study the 
bound state problem by solving the thermodynamic T - matrix for quarkonia 
[5] , or the equivalent Schrodinger- type wave equation where medium effects 
are absorbed in a plasma Hamiltonian (UJ[3]. 

On the other hand one may calculate proper correlators directly from 
lattice QCD and extract from them spectral functions [7j. There is an in- 
triguing disagreement between the Mott temperatures deduced from these 
spectral functions and those of the potential models: From the lattice data 
for quarkonium correlators one has extracted 7j^ l tt ~ 1.9T C while in po- 
tential model calculations Tjj^ ~ 1.2T C . This problem has lead to the 
discussion of the proper thermodynamical function to be used as a poten- 
tial in the Schrodinger equation, see [5] and references therein. 

In this contribution we follow the recently suggested [9J modification of 
the standard one-loop calculation of the Debye mass in thermal Quantum 
Field Theory [10\ [TT] in the framework of the Poyakov-Nambu- Jona-Lasinio 
model, now widely used for a microscopic QCD-motivated description of 
mesons in quark matter [12\ I13j. We then solve the Schrodinger equation 
for charmonium and bottomonium states with the plasma Hamiltonian [3] 
corresponding to the screened Cornell potential [14J and obtain the Mott 
dissociation temperatures of these states. 
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2. Debye-screening in a PNJL quark plasma 

Given the static interaction potential V(q), q 2 = |q| 2 , the statically 
screened potential is given by a resummation of one-particle irreducible 
diagrams ("bubble" resummation = RPA) 

V ac (q) = V(q)/[l + F(0;q)/q 2 ] , (1) 

where the longitudinal polarization function F(0;q) = — IToo(0;q) in the 
finite T case can be calculated within thermal field theory as 

00 f d 3 p 

n 00 (^;q) = g 2 T / 7^33 Tr[7°S'$(iu;„,;p)7 S' <I ,(ia; n - iwjjp- q)] . 

n=— oo J \ I 

(2) 

Here ui = 2ttIT are the bosonic and uo n = (2n + 1)ttT are the fermionic 
Matsubara frequencies of the imaginary-time formalism. The symbol Tr 
stands for traces in color, flavor and Dirac spaces. S§ is the propagator of a 
massless fermion coupled to the homogeneous static gluon background field 
(f3 . Its inverse is given by [121 S3] 

S^ipiWn) =7-p + 7o*w n - A 3 </> 3 , (3) 

where ips is related to the Polyakov loop variable defined by [12] 

$(T) = ^TreCe^^) = 1 (1 + 2cos(/%> 3 )) ■ 

The physics of $(T) is governed by the temperature-dependent Polyakov 
loop potential IA{Q), which is fitted to describe the lattice data for the 
pressure of the pure glue system [12]. After performing the color-, flavor- and 
Dirac traces and making the fermionic Matsubara summation, we obtain in 
the static, long wavelength limit 

n 00 (q) = r d PP 2 ^ = -2g 2 T 2 m = —m 2 D (T) , (4) 

TT Z Jo Op 

where m£)(T) is the Debye mass, the number of degrees of freedom is N c = 3, 
Nf = 2 and /$ (p) is the quark distribution function [13] . For the discussion 
of imaginary parts of the polarization function and their relation to kinetics 
see, e.g., [1 1 j . In comparison to the free fermion case \10\ 111] the coupling 
to the Polyakov loop variable <J>(T) gives rise to a modification of the Debye 
mass, given by the integral 



r/^ 12 f°° j $(1 + 2e- x )e~ x + e~ 3x 

= / dxx — — —. (5) 

~ 9 /o 1 + 3$(1 + e- x )e~ x + e~ 3x 
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The temperature dependence of <J>(T) is taken from Ref. [15] . In the limit 
of deconfinement (<£ = 1), the case of a massless quark gas is obtained 
(1(1) = 1), while for confinement (<£ = 0) one finds that 1(0) = 1/9. 
Taking as the unscreened vacuum potential the one-gluon exchange form 
V(q) = —Ana/q 2 , a = g 2 / (Sir), the Fourier transform of the Debye potential 
results as statically screened potential, V sc (q) = —4-Ka/[q 2 + m 2 D (T)] . 



3. Schrodinger equation and Mott temperatures 

In order to calculate the temperature dependence of the two-particle 
energies E n [(T) for charmonium and bottomonium states in a PNJL quark 
plasma, we solve the Schrodinger equation 

HV\r;T)<t> nl (r;T) = E nl (T)cj> nl (r;T) , (6) 

for the Hamiltonian [3] 

V 2 

HP\r; T) = 2m Q - am D (T) + V(r; T), (7) 

with the screened Cornell potential [TU [3] 

V(r;T) = -%~^(T)r + _ e -mv(T) r] (g) 

where parameters are fitted to the vacuum spectroscopy of heavy quarkonia 
by a s = 0.471 = 3a/4, a = 0.192 GeV 2 and the heavy-quark masses m c = 
1.94 GeV, rrib = 5.1 GeV. Here we use the Debye mass of the previous 
section with the temperature dependence of <£(T) taken from a nonlocal 
PNJL model [15J. Note that the Hamiltonian ([7|) contains a temperature- 
dependent shift of the continuum edge due to the Hartree selfenergies of 
the heavy quarks in the potential ([8]), which results in a definition of the 
dissociation energies as 

E^ S (T) := 2m Q + — - am D - E nl (T) , (9) 
mo 

and of the Mott temperatures as E^f s (T^ ott ) = 0. 

In the upper panels of Fig. [T]we show the temperature dependence of the 
two-particle energies (masses) of the lowest heavy quarkonia states at rest in 
the medium together with the corresponding two-particle continuum edge 
for charmonia (left) and bottomonia (right). The coupling to the Polyakov 
loop potential leads to a suppression of the quark-antiquark excitations 
responsible for the screening of the heavy-quark potential in the present 
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Fig. 1. Upper panels: Temperature dependence of the two-particle energies 
(masses) of the lowest heavy quarkonia states at rest in the medium together with 
the corresponding two-particle continuum edge for charmonia (left) and bottomo- 
nia (right) with (solid lines) and without (dashed lines) coupling to the Polyakov- 
loop potential. Lower panels: Temperature dependence of the binding energies 
Eb = — E^f s , compared to the available thermal energy of medium particles 
Eth = T ■ For discussion, see text. 

model. This entails a stabilization of the bound states (solid lines) relative 
to the case without that coupling (dashed lines). In the lower panels we 
show the corresponding binding energies E-q = —E^f s . Comparison with 
the available thermal energy of medium particles (-Eth = T) shows that 
Mott dissociation by thermal activation [TB] is possible well below the Mott 
temperatures estimated from the solution of the Schrodinger equation. For 
more detailed discussions of the dissociation kinetics of heavy quarkonia, 
see the recent review [3] and references therein. 

4. Conclusions 

We have applied the methods of thermal field theory to estimate the 
effects of Debye screening on heavy quarkonia bound state formation. In 
order to account for residual effects of confining color correlations in the 
deconfined phase, we have used the PNJL model in the evaluation of the 
one-loop polarization function. As expected, a stabilization of bound states 
in the vicinity of the critical temperature for T > T c is obtained. We have 
solved numerically the Schrodinger equation to derive the Mott criterion for 
bound states of the statically screened Cornell potential and obtained Mott 
temperatures in good agreement with previous results from nonrelativistic 
potential models exploiting lattice QCD singlet free energies as potentials 
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in the Schrodinger equation for heavy quarkonia. This agreement with 
previous results (see, e.g., Ref. [8]) includes also higher quarkonia reso- 
nances. For these states the stabilization effect is most pronounced since 
their Mott temperatures lie in the region of temperatures where the sup- 
pression of quark-antiquark excitations due to the Polyakov-loop potential 
is the largest. Due to the lowering of the dissociation energies for heavy 
quarkonia states, their dissociation become possible by thermal activation 
already well before the Mott temperatures are reached. According to the 
present work, even the tightly bound bottomonium ground state T may get 
"boiled away" at temperatures T ~ 250 MeV, well accessible already at 
RHIC. The detailed discussion of the heavy quarkonia dissociation kinetics 
and its relationship to imaginary parts of the heavy-quark potentials [H] is 
beyond the scope of the present contribution. 
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